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We  have  studied  the  formation  and  growth  of  solid-electrolyte  interphase  (SEI)  for  the  case  of  ethylene 
carbonate  (EC),  dimethyl  carbonate  (DMC)  and  mixtures  of  these  electrolytes  using  molecular  dynamics 
simulations.  We  have  considered  SEI  growth  on  both  Li  metal  surfaces  and  using  a  simulation  framework 
that  allows  us  to  vary  the  Li  surface  density  on  the  anode  surface.  Using  our  simulations  we  have  obtained 
the  detailed  structure  and  distribution  of  different  constituents  in  the  SEI  as  a  function  of  the  distance 
from  the  anode  surfaces.  We  find  that  SEI  films  formed  in  the  presence  of  EC  are  rich  in  LhCCU  and  Li20, 
while  LiOCH3  is  the  primary  constituent  of  DMC  films.  We  find  that  dilithium  ethylene  dicarbonate,  LiEDC, 
is  formed  in  the  presence  of  EC  at  low  Li  surface  densities,  but  it  quickly  decomposes  to  inorganic  salts 
during  subsequent  growth  in  Li  rich  environments.  The  surface  films  formed  in  our  simulations  have  a 
multilayer  structure  with  regions  rich  in  inorganic  and  organic  salts  located  near  the  anode  surface  and 
the  electrolyte  interface,  respectively,  in  agreement  with  depth  profiling  experiments.  Our  computed 
formation  potentials  1.0  V  vs.  Li/Li+  is  also  in  excellent  accord  with  experimental  measurements.  We 
have  also  calculated  the  elastic  stiffness  of  the  SEI  films;  we  find  that  they  are  significantly  stiffer  than  Li 
metal,  but  are  somewhat  more  compliant  compared  to  the  graphite  anode. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

As  the  global  energy  demand  increases,  developing  energy  stor¬ 
age  systems  with  higher  energy  densities  is  becoming  more  and 
more  critical.  Rechargeable  Li-ion  batteries  have  been  widely  used 
in  portable  electronics  due  to  their  high  gravimetric  energy  storage 
and  now  are  being  considered  for  critical  applications  such  as  heavy 
automotives  and  medical  devices.  Currently,  the  negative  electrode 
material  used  in  most  Li-ion  batteries  is  graphite,  which  forms 
lithium-graphite  intercalation  compounds  (Li-GICs)  [1,2].  During 
the  first  charge  of  the  Li-ion  battery,  the  electrolyte  decomposes  at 
the  surface  of  the  graphite  anode.  This  forms  a  film  composed  of 
inorganic  and  organic  electrolyte  decomposition  products,  called 
the  solid  electrolyte  interphase  (SEI)  [3].  While  this  layer  imparts 
kinetic  stability  to  the  electrolyte  against  further  reduction  in  suc¬ 
cessive  cycles  and  thereby  ensures  good  cyclability  of  the  electrode, 
it  leads  to  irreversible  loss  of  Li  and  the  capacity  of  the  battery  [4-8]. 
The  composition,  thickness,  morphology,  and  compactness  of  the 
SEI  all  significantly  affect  battery  performance  [9,10].  The  proper¬ 
ties  of  the  SEI  become  even  more  important  during  cycling  at  high 
rates  and  at  deeper  depth  of  discharge.  In  particular,  during  power 
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cycling,  large  stresses  build  up  in  the  electrode  leading  to  crack¬ 
ing  of  the  SEI  films  and  graphite  is  again  exposed  to  the  electrolyte 
[11].  Growth  of  fresh  SEI  on  the  crack  faces  leads  to  further  loss  of 
capacity  [8].  Therefore  understanding  the  formation  mechanisms, 
composition,  structure  and  mechanical  and  transport  properties  of 
the  SEI  is  critical  for  design  of  long-lived,  high-performance  batter¬ 
ies. 

As  discussed  in  detail  in  a  recent  review  article  on  SEI,  the  picture 
of  a  real  SEI  inside  a  Li-ion  battery  has  always  been  a  blur  [9,10]. 
Models  proposed  by  Peled  et  al.  [12],  Aurbach  et  al.  [13,14],  and 
Edstrom  et  al.  [15,16]  suggest  SEI  to  be  a  dense  layer  of  inorganic 
components  close  to  the  anode,  followed  by  a  porous  organic  or 
polymeric  layer  close  to  the  electrolyte  phase.  From  an  experimen¬ 
tal  point  of  view,  since  most  of  the  characterization  studies  of  the 
SEI  are  carried  out  ex-situ ,  it  is  difficult  to  know  with  a  high  degree 
of  certainty  whether  the  structure  of  the  SEI  is  influenced  by  the 
artifacts  of  the  measurement  and  processing  techniques.  There¬ 
fore  detailed  information  of  the  composition  and  morphology  of 
the  SEI  are  currently  under  debate.  For  example,  Leifer  et  al.  very 
recently  suggested  an  entirely  new  family  of  SEI  products,  which 
have  not  been  identified  in  prior  work  using  XPS  and  FT-IR  [17]. 
Given  the  experimental  difficulty  in  characterizing  the  SEI,  a  natu¬ 
ral  question  to  ask  is  whether  computational  models  can  be  used 
to  address  issues  related  to  the  composition  and  structure  of  the 
SEI. 
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Molecular  simulations  (MD)  using  chemically  accurate  meth¬ 
ods  can  provide  information  on  the  reaction  pathways,  rates  and 
byproducts  of  electrolyte  decomposition  reactions  as  well  the  com¬ 
position  and  spatial  distribution  of  the  reduction  products  in  the  SEI 
films.  However,  most  studies  to  date,  that  employ  quantum  simula¬ 
tions,  have  focused  on  modeling  the  former  [18-20].  More  recently 
ab  intio  MD  studies  have  been  used  to  study  the  early  stages  of 
SEI  formation  [21],  but  given  the  computational  cost,  these  cal¬ 
culations  cannot  handle  the  sizes  of  typical  SEI  films  (typically 
lOnm  or  more).  In  this  work,  we  employ  a  reactive  force  field 
potential,  ReaxFF  [22],  to  study  the  structure  of  the  SEI  formed 
when  a  Li  metal  anode  is  “dipped”  in  an  electrolyte  and  in  the 
case  where  Li  atoms  are  introduced  at  the  surface  of  the  simula¬ 
tion  box  and  allowed  to  interact  with  electrolyte  molecules.  In  the 
ReaxFF  that  we  employ  in  this  work,  atomic  connectivity  is  deter¬ 
mined  by  bond  orders  calculated  from  interatomic  distances  that 
are  updated  every  MD  step.  This  allows  for  bonds  to  break  and 
form  during  the  simulation.  In  order  to  account  for  non-bonded 
interactions  such  as  van  der  Waals  and  Coulomb  interactions  for 
a  system  with  changing  connectivity,  these  interactions  are  cal¬ 
culated  between  every  pair  of  atoms,  irrespective  of  connectivity 
[23]. 

There  are  competing  and  parallel  solvent  reduction  processes, 
which  result  in  deposition  of  a  number  of  organic  and  inorganic 
decomposition  products  on  the  surface  of  anode.  To  ensure  that 
the  reactive  potential  we  have  studied  capture  these  fundamen¬ 
tal  processes,  we  first  simulate  several  key  reactions  leading  to 
primary  and  secondary  decomposition  products.  We  find  that  the 
reaction  pathways  and  products  we  obtain  are  generally  in  agree¬ 
ment  with  first  principles  calculations  and  experiments.  Having 
confirmed  the  ability  of  our  approach  to  capture  the  “unit  pro¬ 
cesses”  leading  to  the  formation  of  the  SEI,  we  proceed  to  simulate 
the  growth  of  the  SEI  when  the  electrolyte  is  exposed  to  Li  anode 
or  when  Li  atoms  are  introduced  at  a  putative  anode  surface.  Using 
these  simulations  we  obtain  the  detailed  structure  of  the  SEI  includ¬ 
ing  its  composition  as  a  function  of  the  distance  from  the  anode 
surface.  Experiments  have  shown  that  the  SEI  compositions  and 
contents  are  a  direct  consequence  of  the  electrolyte  composition. 
We  have  therefore  carried  out  these  simulations  with  Ethylene 
Carbonate  (EC),  Dimethyl  Carbonate  (DMC)  and  mixtures  of  both. 
We  find  that  lithium  alkoxides  appear  in  the  SEI  formed  in  DMC 
electrolytes,  while  Li20  and  Li2C03  appear  in  the  SEI  formed  in  EC 
(Fig.  5).  Both  solvents  react  to  form  several  products,  both  R0C02Li 
and  ROLi  species.  Further  reduction  of  R0C02Li  leads  to  the  for¬ 
mation  of  surface  Li2C03.  Consequently,  the  surface  films  formed 
in  our  simulations  have  a  multilayer  structure  as  confirmed  by 
depth  profiling  and  XPS  [24],  and  are  laterally  non-uniform,  as 
confirmed  by  in  situ  AFM  [25].  Under  the  same  conditions,  the 
thickness  of  SEI  formed  with  EC  is  more  than  DMC,  indicating 
that  the  EC  is  more  reactive,  also  in  accord  with  observations 
[26,27].  Our  simulations  also  allow  us  to  compute  the  formation 
potential  of  the  SEI  for  different  electrolytes.  Finally,  a  good  SEI 
should  be  a  compact  layer  adhering  well  to  the  anode.  It  should 
be  elastic  and  flexible  to  accommodate  non-uniform  electrochem¬ 
ical  behavior  and  active  material  breathing.  At  present,  there  is  no 
information  on  the  mechanical  properties  of  the  SEI  and  how  these 
properties  are  related  to  the  composition  of  the  electrolyte.  We 
have  used  our  simulations  to  study  the  elastic  properties  of  the 
SEI. 


2.  Simulation  methods 

MD  simulations  were  performed  using  ReaxFF,  which  is  a 
general  bond-order-dependent  potential  that  provides  accurate 
descriptions  of  bond  breaking  and  bond  formation  [22].  The  main 


Fig.  1.  (a)  Schematic  of  a  Li  metal  electrode  dipped  in  an  electrolyte,  (b)  Initial 
configuration  of  the  cell  used  in  our  simulations. 


difference  between  traditional  nonreactive  force  fields  and  the 
one  that  we  employ  in  this  work  is  that  the  atomic  connectiv¬ 
ity  is  determined  by  bond  orders  calculated  from  interatomic 
distances  that  are  updated  every  MD  step.  This  allows  for 
bonds  to  break  and  form  during  the  simulation.  In  order  to 
account  for  non-bonded  interactions  such  as  van  der  Waals  and 
Coulomb  interactions  for  a  system  with  changing  connectivity, 
these  interactions  are  calculated  between  every  pair  of  atoms, 
irrespective  of  connectivity,  and  any  excessive  close-range  non- 
bonded  interactions  are  avoided  by  the  inclusion  of  a  shielding 
term  [23].  The  force  field  used  in  the  current  work  has  been 
trained  against  the  quantum-mechanical  data  for  Li-C-O-H  sys¬ 
tems  described  in  Han  et  al.  [28]  and  was  re-optimized  recently. 
As  we  show  below,  this  potential  can  provide  the  details  of  the 
variations  of  the  composition  of  the  SEI  over  nanometer  length 
scales. 

To  model  the  formation  of  SEI,  we  have  considered  two  types  of 
simulations.  In  the  first  case,  we  model  the  fundamental  electrolyte 
decomposition  reactions  that  occur  when  Li  metal  is  “dipped”  in  an 
electrolyte  [12,29-31]  as  shown  in  Fig.  1(a).  Here,  as  the  Li  atoms 
interact  with  the  electrolyte,  they  transfer  electrons  to  the  elec¬ 
trolyte  molecules  forming  organic  and  inorganic  decomposition 
products.  As  we  show  below,  the  advantage  of  these  simulations  is 
that  basic  decomposition  reactions  can  be  studied  without  the  need 
for  “charging”  the  anode.  The  electrolyte  at  the  equilibrium  density 
of  liquid  phase  is  filled  above  single  crystal  bcc  Li  with  dimensions 
(8  x  16  x  20  <3q).  The  systems  we  consider  are  composed  of  5120  Li 
atoms  and  about  800  electrolyte  molecules.  The  size  of  the  simula¬ 
tion  box  is  approximately  3  nm  x  6  nm  x  14  nm.  Periodic  boundary 
conditions  are  enforced  along  the  x  and  y  directions  (Fig.  1(b)). 
In  the  z  direction,  a  flat  wall,  which  interacts  with  the  atoms  by 
generating  a  force  in  a  direction  perpendicular  to  the  wall  to  pre¬ 
vent  them  from  escaping  during  MD  simulations,  is  introduced.  In 
our  simulations,  charge  transfer  is  performed  by  the  charge  equi¬ 
libration  (QEq)  method  [32]  at  every  MD  step.  The  temperature  of 
the  system  is  controlled  by  rescaling  atomic  velocities  every  10 
MD  steps  (each  MD  time  step  At=0.2fs).  The  system  is  heated 
up  to  the  desired  temperature  during  the  first  5  x  104  MD  steps 
(=  10  ps)  and  kept  at  constant  temperature  for  40  ps.  Based  on  the 
bond  information  computed  by  the  ReaxFF  potential,  a  home-built 
program  was  developed  to  that  identify  the  constituents  of  the  SEI 
based  on  the  number  of  bonds  and  charge  state  of  each  atom  in  a 
molecule. 

In  the  second  set  of  calculations,  we  simulate  SEI  formation  dur¬ 
ing  the  charging  process  as  Li+  ions  from  the  electrolyte  combine 
with  the  electrons  from  the  anode  and  the  electrolyte  molecules  in 
the  vicinity  of  the  anode  surface  (typically  less  than  a  nanometer, 
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Fig.  2.  Schematic  depiction  of  the  SEI  formation  during  initial  charging  cycles. 


which  is  the  distance  over  which  the  electron  can  tunnel  from  the 
anode  [33]).  In  our  molecular  simulations,  it  is  difficult  to  confine 
the  electrons  to  the  anode  and  Li+  ions  to  the  electrolyte.  How¬ 
ever,  since  the  electrons  and  Li+  ions  combine  with  the  electrolyte 
to  form  the  SEI  only  in  the  immediate  vicinity  of  the  anode,  this 
difficulty  can  be  circumvented  by  introducing  Li  atoms  (or  equiva¬ 
lently  Li+  +  e-)  within  a  distance  of  1  nm  from  the  anode  surface 
(Fig.  2).  The  atoms  can  combine  with  the  electrolyte  molecules 
to  form  the  SEI  products,  thus  mimicking  the  charging  process.  In 
our  simulations,  we  introduce  a  certain  number  of  Li  atoms  at  the 
anode  surface  and  allow  them  to  react  with  electrolyte  molecules. 
Li  atoms  are  introduced  every  1 .0  ps  at  randomly  selected  locations 
on  the  bottom  surface  of  the  simulation  cell  ( x-y  plane),  whose  size 
is  3.3  nm  x  3.4  nm  x  4.5  nm  in  x,y,  and  z  directions,  and  the  temper¬ 
ature  of  the  system  is  kept  at  300  K  by  rescaling  the  temperature 
every  10MD  steps.  This  rate  of  introducing  the  Li  atoms  is  suffi¬ 
cient  to  observe  the  initial  stages  of  the  SEI  growth  because  the 
major  reactions  such  as  bond  breaking  and  reformation  between  Li 
ions  and  electrolytes  occur  within  a  few  ps  [21  ].  As  in  the  case  of  the 
first  set  of  calculations,  periodic  boundary  conditions  are  enforced 
in  x  and  y  directions.  We  note  that  the  typical  concentration  or 
areal  density  of  Li+  ions  at  the  surface  of  C  anode  is  3.63/nm2  for 
zigzag  edge  or  6.28/nm2  for  armchair  edge,  which  is  significantly 
smaller  that  the  concentration  of  Li  atoms  on  the  surface  of  Li  metal 
(9.3/nm2).  As  we  discuss  below,  the  concentration  of  Li  atoms  has 
a  significant  influence  over  the  type  of  the  SEI  products  that  are 
formed. 


4-.*  4“  -  V* 


Fig.  3.  Basic  decomposition  reactions  for  the  case  of  (a)  ethylene  carbonate  (EC), 
and  (b)  dimethyl  carbonate  (DMC)  interacting  with  Li.  The  dashed  lines  indicate  the 
bonds  that  are  broken  as  a  result  of  the  reactions. 


3.  Results  and  discussion 

3.1.  Mechanisms  for  the  formation  of  SEI  components 

The  structure  of  the  electrolyte  molecules  (EC  and  DMC)  relaxed 
using  the  ReaxFF  potential  is  in  the  first  frame  shown  in  Fig.  3. 
Comparing  the  computed  structure  with  density  functional  the¬ 
ory  (DFT)  calculations  and  experiments  (Table  1),  we  find  that 
the  potential  captures  all  the  bond  lengths  and  angles  very  well. 
Next,  we  consider  the  key  decomposition  reactions  that  involve 
interaction  of  the  electrolyte  with  Li  atoms.  To  study  these  reac¬ 
tions,  we  introduce  Li  atoms  randomly  at  distances  >5  A  from  the 
centroids  of  the  electrolyte  molecules  and  allow  the  molecules 
to  interact  at  temperatures  close  to  500 1<  for  2  ps.  We  tried  sev¬ 
eral  sets  of  initial  positions,  but  the  results  predominantly  fall  into 
the  two  categories  shown  in  Fig.  3.  As  in  the  case  of  DFT  calcula¬ 
tions  [21],  the  reactions  occur  in  most  cases  within  1.0  ps  or  less. 
Fig.  3(a)  shows  the  primary  reaction  pathways  for  a  Li  atom  inter¬ 
acting  with  an  EC  molecule.  The  reaction  results  in  the  formation 
of  a  Li0C02C2H4  radical  by  means  of  breaking  the  CE-Oi  bond, 
in  agreement  with  previous  DFT  studies  [18].  After  the  cleavage 
of  this  bond,  the  charge  of  the  atoms,  CE  and  Oi  changed  from 
+0.01 6e  and  -0.47e  to  -0.23e  and  -0.61  e,  respectively.  As  we  dis¬ 
cuss  below,  Li0C02C2H4  can  further  react  to  from  molecules  such 
as  dilithium  ethylene  dicarbonate  ([Li0C02CH2]2:LiEDC),  lithium 
carbonate  (Li2C03)  and  Li20.  Two  primary  pathways  for  the  inter¬ 
action  of  Li  atoms  with  DMC  are  shown  in  Fig.  3(b).  In  one  case, 
lithium  methoxide  (LiOCH3)  is  formed  as  a  result  of  the  cleavage  of 
the  Cc-Oi  bond,  while  the  other  reaction  leads  to  the  formation  of 
lithium  methyl  carbonate  (Li0C02CH3)  as  a  result  of  the  cleavage 
of  the  CM-Oi  bond.  The  resulting  molecule,  OCOCH3  is  the  precur¬ 
sor  for  the  formation  of  CO  and  LiOCH3  as  discussed  below.  These 
results  for  the  decomposition  of  DMC  are  in  good  agreement  with 
the  reaction  suggested  on  the  basis  of  experimental  measurements 
[36,37]. 

Next  we  focus  on  secondary  decomposition  reactions  that 
involve  further  interactions  of  the  primary  decomposition  products 
in  Fig.  3  with  each  other  or  with  Li.  The  reaction  pathway  for  the 
formation  of  LiEDC,  Fig.  4(a),  and  a  hydrocarbon  gas  molecule,  C2H4 
as  a  result  of  the  interaction  between  two  Li0C02C2H4  radicals  is 
shown  in  Fig.  3(a).  Generally,  LiEDC  is  considered  to  be  a  primary 
component  of  the  SEI  in  EC  based  electrolytes  [38,39].  However, 
we  find  that  this  molecule  can  easily  breakdown  in  the  presence  of 
Li  atoms.  As  shown  in  Fig.  4(b),  Li  atoms  directly  attack  the  Oi-CE 
bond  in  LiEDC  and  dissociate  it  into  lithium  carbonate  (Li2C03 )  and 
Li0C02C2H4.  The  latter  radical  can  further  interact  with  Li  atoms  to 
eventually  form  two  lithium  carbonate  molecules  and  an  ethylene 
gas  molecule.  Lithium  carbonate  can  further  decompose  (Fig.  4(c)) 
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Table  1 

Characteristics  of  ethylene  carbonate  (upper  values)  and  dimethyl  carbonate  (lower  values).  The  units  of  the  length,  charge  and  angle  are  A,  e  and°,  respectively. 


This  work 

Brown  [34] 

Wang  et  al.  [18] 

Tasaki  et  al.  [19] 

Xing  et  al.  [35] 

02=Cc 

1.152 

1.15 

1.187 

- 

1.19 

1.159 

- 

- 

- 

1.21 

O 

1 

n 

n 

1.495 

1.33 

1.353 

- 

1.36 

1.460 

- 

- 

- 

1.34 

Oi-Ce 

1.506 

1.40 

1.426 

- 

1.44 

Oi-Cm 

1.458 

- 

- 

- 

1.44 

qOi 

-0.469 

- 

- 

-0.36 

-0.55 

-0.489 

- 

- 

-0.39 

- 

q02 

-0.533 

- 

- 

-0.50 

-0.56 

-0.616 

- 

- 

-0.55 

- 

qCc 

0.955 

- 

- 

0.80 

1.01 

0.961 

- 

- 

0.92 

- 

qCE 

0.016 

- 

- 

0.01 

-0.05 

qC  M 

-0.061 

- 

- 

-0.15 

- 
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- 
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- 

- 

- 
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in  the  presence  of  Li  atoms  to  form  Li20  and  a  CO  molecule.  Pres¬ 
ence  of  CO  has  been  detected  in  recent  experimental  studies  [40]. 
Fig.  4(d)  shows  the  mechanism  for  the  formation  of  CO  and  a  LiOCH3 
molecule  from  the  OCOCH3  radical,  which  is  one  of  the  primary 
structures  that  results  from  the  reaction  between  a  DMC  and  a  Li 
(Fig.  3(b)).  Finally,  Fig.  4(e)  shows  the  formation  mechanism  of  a 
Li2C03  molecule  from  the  Li0C02CFI3  radical.  Having  verified  that 
our  MD  simulations  predict  primary  and  secondary  decomposi¬ 
tion  reactions  in  agreement  with  available  DFT  calculations  [18] 
and  pathways  deduced  from  experiments  [36],  we  now  proceed  to 
study  the  structure  of  the  SEI  when  all  of  these  reactions  can  simul¬ 
taneously  occur.  Note  that  given  the  computational  cost,  quantum 
mechanical  approaches  cannot  be  used  to  determine  the  structure 
of  the  ~10nm  thick  SEI  films,  but  as  we  show  next,  the  molecu¬ 
lar  simulations  we  employ  can  predict  the  spatial  variations  in  the 
composition  of  the  constituents  of  SEI,  their  formation  potentials 
and  elastic  properties. 


b  Jt*  *  ^ 


Fig.  4.  Secondary  decomposition  reactions  for  the  key  SEI  components,  (a)  Forma¬ 
tion  of  LiEDC  and  a  C2H4  molecule  from  two  UOCO2C2H4  radicals,  (b)  Two  Li2C03 
and  a  C2H4  formed  from  a  LiEDC  molecule  interacting  with  additional  Li  atoms,  (c) 
Two  Li20  and  a  CO  molecule  formed  from  Li2C03.  (d)  A  LiOCH3  and  a  CO  molecule 
formed  from  OCOCH3.  The  dashed  lines  are  the  bonds  that  are  broken  as  a  result  of 
the  reactions.  Dark  gray,  white,  red  and  big  light  gray  atoms  denote  C,  H,  0  and  Li, 
respectively. 


3.2.  Effects  of  electrolytes  on  the  evolution  of  the  SEI 

The  atomic  configurations  after  40  ps  NVT-MD  simulation  at 
300  K  and  500  K  for  the  electrolyte-Li  metal  system  are  shown  in 
Fig.  5.  The  constituent  atoms  are  color-coded  according  to  their 
charge,  which  allows  us  to  readily  distinguish  the  SEI  region  from 
the  electrode;  Li  atoms  are  neutral  in  the  electrode  and  are  in  a  posi¬ 
tive  charge  state  (>+0.1  e)  once  they  interact  with  the  electrolyte  and 
become  part  of  the  SEI.  Our  simulations  (Fig.  5(b))  show  that  the  SEI 
films  grow  faster  in  the  case  of  EC  compared  to  DMC,  with  EC  +  DMC 
mixtures  falling  in  between.  As  shown  in  Fig.  5(c),  the  density  of  Li 
consumed  in  the  SEI  tends  to  an  asymptotic  limit,  which  depends 
on  the  type  of  electrolyte  but  weakly  on  the  temperature.  The  satu¬ 
ration  values  lie  between  0.1  and  0.13  g/cm3,  which  is  much  lower 
than  the  density  of  Li  metal  (0.53  g/cm3).  The  density  of  Li  incor¬ 
porated  in  the  case  of  EC  is  more  than  that  amount  consumed  in 
DMC  and  the  case  of  EC  +  DMC  mixtures  lies  in  between.  This  result 
agrees  with  the  experimental  observations  that  show  that  EC  is 
more  reactive  compared  to  DMC  [26,27]. 

In  the  case  of  Li  metal,  the  constituents  of  the  SEI  components  are 
highly  dependent  on  the  type  of  the  electrolyte  as  shown  in  Fig.  6. 
The  SEI  in  the  case  of  EC  is  composed  largely  of  Li20  and  Li2C03 
along  with  trapped  gas  molecules,  C2H4  and  CO.  We  find  that  while 
LiEDC  is  formed  in  the  very  early  stages  of  SEI  growth,  it  decomposes 
into  other  byproducts  by  interacting  with  Li  atoms  of  the  metal 
anode.  Since  LiEDC  is  considered  to  be  a  major  component  of  SEI  in 
graphite  anodes,  we  hypothesize  that  it  would  not  decompose  in 
the  Li-poor  environment  found  at  the  surface  of  graphite  compared 
to  the  case  of  Li  anode.  To  test  this  hypothesis,  we  carried  out  SEI 
growth  simulations  using  the  approach  sketched  in  Fig.  2,  where 
Li  atoms  are  introduced  at  the  putative  anode  surface.  Indeed,  as 
shown  in  Fig.  7  we  find  copious  formation  of  LiEDC  and  Li2C03  at 
early  stages  of  SEI  growth.  As  the  growth  of  the  SEI  progresses,  the 
quantity  of  Li2C03  keeps  increasing,  but  the  number  of  the  LiEDC  in 
this  case  peaks  around  5.0atoms/nm2  and  subsequently  saturates. 
The  amount  of  Li20  in  the  SEI  increases  at  the  expense  of  LiEDC.  In 
the  case  of  Li  metal  anode,  given  the  high  concentration  of  Li  atoms, 
almost  all  of  the  LiEDC  that  is  formed  is  instantaneously  consumed. 
Taken  together,  the  results  in  Figs.  6  and  7  show  that  Li  concentra¬ 
tion  plays  an  important  role  in  determining  the  composition  of  SEI 
in  the  case  of  EC.  For  the  case  of  high  Li  concentration  as  in  Li  metal, 
the  bond  between  carbonate  (C03)  and  ethylene  (C2H4)  can  easily 
be  broken  by  Li  atoms  leading  to  formation  of  hydrocarbon  gases. 
The  carbonate  group  that  is  left  behind,  reacts  with  additional  Li 
ions  and  consequently  dissociates  into  the  different  types  of  Li  salt 
such  as  Li2C03  and  Li20  as  shown  in  Fig.  4(c).  Very  recent  exper- 
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Fig.  5.  (a)  Atomic  configuration,  (b)  evolution  of  the  thickness  of  the  SEI  and  (c)  Li  density  in  the  SEI  formed  between  Li  metal  and  electrolytes  (EC,  EC  +  DMC,  and  DMC)  at 
different  temperatures  (300  K  and  500 1<).  Snapshots  in  (a)  are  obtained  after  40  ps  of  NVT  MD  simulations.  The  atoms  are  color-coded  according  to  their  charge  state  (green 
is  neutral). 


Fig.  6.  Distribution  of  the  SEI  components  for  different  electrolytes  (left).  Atomic  configurations  from  MD  simulations;  the  components  of  the  SEI  are  identified  (right). 
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Fig.  7.  The  evolutions  of  EC  based  SEI  components  with  increasing  Li  density.  Inset 
is  a  snapshot  of  the  atomic  configurations  when  the  Li  (blue  atoms)  surface  density 
is  20/nm2. 


imental  work  on  the  composition  of  SEI  [41  ]  on  electrodeposited 
Li  metal  indeed  shows  that  Li2C03  and  Li20  are  the  primary  com¬ 
ponents  of  the  SEI  in  agreement  with  our  calculations.  While  this 
work  found  the  presence  of  LiOH  and  Li3N  in  SEI  films  (due  to  expo¬ 
sure  to  air),  LiEDC  was  not  detected  in  the  SEI  films  confirming  the 
predictions  of  our  work. 

The  SEI  in  the  case  of  DMC  is  composed  of  Li20,  Li0C02CH3, 
Li2C03,  and  LiCH3,  with  LiOCH3  being  the  dominant  constituent 
[38,39].  If  the  dissociation  follows  the  first  of  the  reaction  path¬ 


ways  shown  in  Fig.  3(b),  two  LiOCH3  molecules  can  be  formed  along 
with  a  CO  molecule.  However,  if  a  Li  atom  reacts  with  the  carbon¬ 
ate  part  of  the  DMC  molecule  (which  corresponds  to  the  cleavage 
of  the  Cm-Oi  bond  as  shown  in  Fig.  3(b)),  the  types  of  the  byprod¬ 
ucts  that  are  formed  are  determined  by  the  number  of  Li  atoms 
that  replace  the  methyl  radicals.  Li2C03  is  created  when  two  Li 
atoms  replace  two  methyl  radicals  (Fig.  4(e)),  while  Li0C02CH3  is 
formed  when  only  methyl  radical  is  replaced  (Fig.  4(d)).  As  in  the 
case  of  EC,  where  C2H4  molecules  are  formed  during  decomposi¬ 
tion,  we  can  expect  the  formation  of  ethane  (C2H6)  from  DMC,  but 
we  did  not  detect  it  to  any  significant  extent.  This  is  because  of 
the  small  probability  for  the  completion  of  the  reaction  between 
the  two  methyl  radicals,  *CH3  +  #CH3  C2Hg,  which  is  needed  for 
the  formation  of  ethane.  Instead,  LiCH3  is  formed  when  Li  atoms 
combine  with  the  methyl  radical  [42].  In  DMC,  CO  turns  out  to  be 
the  major  gas  molecule  produced  (which  constitutes  31.6%  of  the 
total  amount  of  the  SEI  components),  in  good  agreement  with  the 
previous  experimental  work  [43].  When  Li  reacts  with  the  mixture 
electrolyte  of  EC  and  DMC,  the  type  of  SEI  components  follow  the 
order  Li20  >  LiOCH3  >  Li2C03  >  Li0C02CH3  >  LiCH3.  The  reason  why 
Li20  is  the  most  abundant  component  is  because  the  secondary 
reactions  of  Li  atoms  with  both  EC  and  DMC  decomposition  prod¬ 
ucts  lead  to  the  formation  of  Li20  molecules. 

Next,  we  consider  the  spatial  distribution  of  the  SEI  components 
in  Fig.  8.  The  surface  films  formed  in  our  simulations  have  a  mul¬ 
tilayer  structure  in  agreement  with  depth  profiling  experiments. 
We  find  that  the  gases  formed  as  a  result  of  electrolyte  decompo¬ 
sition  reactions  are  trapped  near  the  Li  metal  side.  Given  the  short 
time  scale  of  our  simulations,  these  gases  are  not  able  to  diffuse 
through  the  SEI;  in  real-life,  these  molecules  can  escape  from  the 
SEI-electrode  interface  by  diffusion.  Inorganic  Li  salts  such  as  Li2C03 
and  Li20  tend  to  be  closer  to  the  Li  metal  electrode,  while  organic  Li 
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Fig.  8.  Spatial  distribution  of  the  SEI  components  for  different  electrolytes. 
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Fig.  9.  The  formation  potential  for  different  electrolytes  plotted  as  a  function  of  the  thickness  of  SEI  films. 


salts  such  as  LiOCH3  and  Li0C02CH3  are  found  at  the  SEI-electrolyte 
interface  in  agreement  with  experiments  [  1 2  ].  The  reason  why  inor¬ 
ganic  salts  with  high  Li  content  compared  to  organic  salts  are  found 
near  the  electrode  is  that  these  salts  and  hydrocarbon  gases  result 
from  the  decomposition  of  the  primary  SEI  products  as  they  com¬ 
bine  with  additional  Li  atoms.  High  Li  content  at  the  surface  of  the 
anode  allows  these  reactions  to  proceed  at  the  higher  rate  com¬ 
pared  to  the  electrolyte-SEI  interface,  which  is  dominant  in  primary 
decomposition  products. 

3.3.  Effects  of  electrolytes  on  the  SEI  properties 

Our  simulations  show  that  Li  reacts  spontaneously  with  the 
electrolyte  to  produce  decomposition  products  which  are  more 
thermodynamically  stable.  The  difference  in  the  energy  of  the  reac¬ 
tants  and  the  products  gives  the  so-called  “SEI  formation  potential”. 
In  order  to  quantify  the  formation  potential,  we  use  the  potential 
energy  computed  as  a  function  of  the  SEI  thickness  in  the  equation 
Ef(d)  =  ^^^ctroiyte),  where  NLi+  is  the  number  of  Li+  ions  that 

have  reacted  to  form  the  SEI  and  E(d)  is  the  potential  energy  of  the 
SEI  as  a  function  of  its  thickness  d,  ELi  denotes  the  energy  of  Li  atom 
in  metal  and  Eeiectroiyte  is  the  potential  energy  of  the  electrolyte 
prior  to  decomposition. 

For  the  electrolyte-Li-metal  system,  Fig.  9  shows  the  SEI  forma¬ 
tion  potential,  Ef,  as  a  function  of  the  SEI  thickness  d  for  EC,  DMC 
and  EC  +  DMC  for  different  temperatures.  In  all  cases,  the  forma¬ 
tion  potential  tends  to  an  asymptotic  limit,  which  depends  weakly 
on  temperature.  The  saturation  values  lie  around  1.0  V  vs.  Li/Li+, 
in  good  agreement  with  the  experimental  measurements  that  are 
generally  in  the  range  0.8-1. 0  V  vs.  Li/Li+  [44-47].  The  potential  for 
DMC,  EC  and  EC  +  DMC  are  0.9  V  vs.  Li/Li+,  1 .1  V  vs.  Li/Li+  and  1 .0  V 
vs.  Li/Li+,  respectively. 

A  good  SEI  should  be  elastic  and  flexible  to  accommodate  non- 
uniform  electrochemical  behavior  and  active  material  breathing. 
To  estimate  the  mechanical  compliance  of  the  SEI,  we  now  turn 
to  the  calculation  of  the  elastic  constants  of  the  SEI.  To  achieve 
this  task,  we  first  extract  the  SEI  region  from  the  simulation  box 
excluding  both  Li  metal  and  electrolyte  region  as  shown  schemati¬ 
cally  in  Fig.  10.  The  extracted  sample  is  then  equilibrated  using  NPT 
molecular  dynamics  at  300  K  and  1  atm  for  40  ps.  The  elastic  con¬ 
stants  ( Cjj )  are  obtained  by  calculating  the  change  in  stress  as  small 
strains  ranging  from  0.1  to  1.0%  are  applied  to  the  sample  along  the 
coordinate  axes.  From  the  knowledge  of  the  elastic  constants  (Cy ), 
bulk  modulus  ( B )  and  shear  modulus  (G)  are  then  calculated  using 
the  equation 

Cn  +  C22  +  c33  +  2(Ci2  +  Ci3  +  C23) 


G  _  Cn  +  C22  +  C33  -  Ci 2  -  Ci3  -  C23  C44  +  C55  +  C66 
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Fig.  10.  The  schematic  of  the  method  used  for  computing  the  elastic  properties  of 
the  SEI. 


Table  2 

Mechanical  properties  of  SEI  films. 


SEIec 

SEIec+dmc 

SEIdmc 

Li  metal 

Graphite  [48] 

E  (GPa) 

18.54 

15.99 

14.96 

4.9 

32.47 

B  (GPa) 

14.64 

10.24 

11.47 

11 

29.29 

The  Young’s  modulus  is  then  computed  by  assuming  isotropic 
linear  elastic  stress-strain  relations.  Table  2  shows  the  elastic  prop¬ 
erties  computed  in  the  three  cases.  According  to  our  results,  the 
Young’s  modulus  is  somewhat  larger  for  the  SEI  formed  in  the  case 
of  EC  compared  to  the  case  of  DMC.  Note  also  that  the  SEI  films  are 
more  compliant  to  typical  anode  materials  (graphite  E  =  32.47  GPa), 
but  is  significantly  stiffer  than  Li.  The  elastic  properties  we  have 
computed  here  can  be  used  in  continuum  models  for  stress  gener¬ 
ation  in  Li-ion  battery  anodes. 

4.  Summary  and  conclusions 

In  summary,  we  have  used  molecular  dynamics  simulations 
using  the  ReaxFF  potential  to  simulate  the  formation  of  SEI  on  Li 
metal  anodes  and  by  introducing  Li  atoms  randomly  at  the  surface 
of  our  simulation  cell.  The  latter  set  of  simulations  is  used  to  assess 
the  role  of  Li  density  on  the  formation  of  SEI  films.  While  quantum 
simulations  can  accurately  predict  the  kinetics  of  individual  reac¬ 
tions,  they  cannot  be  used  to  study  competing  reactions  leading 
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to  the  formation  of  complex  SEI  films  typically  several  nanome¬ 
ters  thick.  In  our  MD  simulations,  we  first  considered  the  primary 
and  secondary  decomposition  reactions  that  involve  Li  and  the 
electrolyte  molecules  and  showed  that  our  potential  correctly  cap¬ 
tures  the  reaction  pathways  and  byproducts.  We  then  simulated 
the  formation  of  SEI  allowing  electrolytes  composed  of  EC,  DMC 
and  EC  +  DMC  with  Li  atoms  from  at  the  anode  surface.  Using  our 
simulations,  we  studied  the  distribution  of  organic  and  inorganic 
salts  as  a  function  of  the  distance  from  the  anode  surface.  We  find 
that  Li2C03  and  Li20  are  the  primary  components  of  the  SEI  in 
Li  metal  anodes  when  EC  is  used  as  an  electrolyte.  LiEDC  that  is 
formed  in  the  early  stages  of  SEI  formation  decomposed  into  inor¬ 
ganic  salts  and  hydrocarbon  molecules  in  a  Li  rich  environment.  In 
the  case  of  DMC,  we  find  LiOCH3  to  be  the  dominant  component 
of  the  SEI,  in  excellent  agreement  with  experiments  [37].  Our  sim¬ 
ulations  show  that  inorganic  salts  are  found  closer  to  the  anode 
surface  while  the  region  near  the  electrolyte-SEI  interface  is  rich 
in  organic  salts.  Using  our  simulations  we  find  the  SEI  formation 
potential  to  be  in  the  range  0.9-1. OV  vs.  Li/Li+,  also  in  accord  with 
experiments.  Finally,  we  have  computed  the  elastic  properties  of 
the  SEI;  SEI  is  found  to  be  significantly  stiffer  than  Li  metal,  but  is 
somewhat  softer  than  the  graphite  anode.  Cracking  of  SEI  during 
power  cycling  is  one  of  the  causes  for  the  loss  of  capacity  of  Li-ion 
batteries.  Work  is  currently  underway  to  include  the  properties 
computed  here  in  continuum  models  for  progressive  crack  propa¬ 
gation  in  electrode  particles  during  charge  and  discharge.  Finally, 
while  Li  salts  are  known  to  play  an  important  role  in  the  formation 
of  SEI,  we  have  only  studied  the  decomposition  of  the  bare  elec¬ 
trolyte.  Development  of  new  reactive  potentials  for  P  and  F  will 
allow  one  to  extend  the  present  work  to  account  for  the  presence 
of  salts. 
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